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We study the temperature dependence of the pseudoscalar meson properties in a relativistic bound- 
state approach exhibiting the chiral behavior mandated by QCD. Concretely, we adopt the Dyson- 
Scliwinger approach with a rank-2 separable model interaction. After extending the model to the 
strange sector and fixing its parameters at zero temperature, T = 0, we study the T-dependence 
of the masses and decay constants of all ground-state mesons in the pseudoscalar nonet. Of chief 
interest are i] and rj' . The influence of the QCD axial anomaly on them is successfully obtained 
through the Witten-Veneziano relation at T = 0. The same approach is then extended to T > 0, 
using lattice QCD results for the topological susceptibility. The most conspicuous finding is an 
increase of the rj' mass around the chiral restoration temperature Tch, which would suggest a 
suppression of r]' production in relativistic heavy-ion collisions. The increase of the rj' mass may 
also indicate that the extension of the Witten-Veneziano relation to finite temperatures becomes 
unreliable around and above Tch- Possibilities of an improved treatment are discussed. 

PACS numbers: ll.lO.St, ll.lO.Wx, 12.38.-t, 12.38.Aw, 14.40.Aq 
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I. INTRODUCTION AND PRELIMINARIES 

The experiments at RHIC advanced the empirical 
knowledge about relativistic hot QCD matter dramati- 
cally [H, Q but its properties arc so intricate, that the un- 
derstanding of the state of matter that has been formed is 
still to be developed . Prior to its observation at RHIC, 
the hot QCD matter, usually called quark-gluon plasma 
(QGP), was pictured as a perturbatively interacting gas 
of deconfined quarks and gluons. It was expected to be 
reached above the critical temperature Tc ~ 170 MeV 
0, 0, H, Q where the rapid increase of the effective num- 
ber of active degrees of freedom [e(T)/T^] was found in 
lattice QCD simulations 0, [ll- Nowadays it is clear that 
such a state will be reached only at significantly higher 
temperatures than those accessible today. For tempera- 
tures a few times higher than the critical temperature Tc, 
the interactions and correlations in the hot QCD matter 
are still strong (e.g., see Refs. 1,13), the recent termi- 
nology for it being strongly coupled QGP (sQGP) p^ . 

Critical assessments of hot QCD physics, e.g. Ref. [2], 
stress the present absence of, and the need for, a direct, 
compelling "smoking gun" signal for production of a new 
form of matter in RHIC collisions. The most compelling 
"smoking gun" would be a restoration of the symmetries 
of the QCD Lagrangian in hot, dense matter, notably the 
[SUa(3) flavor] chiral symmetry and the Ua(1) symme- 
try, which are broken in the vacuum. 

Among the issues pointed out as important was also 
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the need to clarify the role of quark-antiquark {qq) bound 
states continuing existence above the critical temper- 
ature Tc 01 ■ Namely, evidence is accumulating that 
strong correlations in the form of quark-antiquark {qq) 
bound states and resonances still exist 0, [Tl| in the 
sQGP above the critical temperature. While in the old 
paradigm even deeply bound charmonium (cc) states 
such as J/'^ and r]c were expected to dissociate at 
T K, Tc, the behavior of their spectral functions as ex- 
tracted by the maximum entropy method (l^ from lat- 
tice QCD simulations of mesonic correlators now indi- 
cates that they should persist till around 2Tc [H, [l3l or 
even above [15|. There are similar indications for light- 
quark mesonic bound states from lattice QCD [3l and 
from other methods 0, ■ This motivates their study 
in the framework of relativistic bound-state equations, 
i.e., the Poincare-covariant Dyson-Schwinger (DS) ap- 
proach to quark- hadron physics [l^, [I^, [l^, [IH, which 
is a continuous approach complementary to lattice QCD 
studies. The DS approach is especially valuable in the 
light-quark sector, where chiral symmetry is essential. 
Namely, the DS approach is unique in that it can incor- 
porate the correct chiral behavior of QCD, unlike other 
bound-state approaches. This is because the crucial low- 
energy QCD phenomenon of dynamical chiral symmetry 
breaking (DChSB) is well-understood and under control 
in DS approach [3j [l^, 113, HH, at least in the consis- 
tent rainbow-ladder approximation (RLA), where kernel 
Ansdtze of the form 
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(where e, f, g,h schematically represent spinor, color and 
flavor indices) are used for the interactions between 
quarks in both Eq. ([2]), the gap equation = S^^ — E 
for the full, dressed quark propagator S {So is the free, 
bare one) and the Bethe-Salpcter (BS) equation ^ for 
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the meson quark-antiquark bound state A4. That is, in 
RLA, the full propagator of the quark of the flavor q 
characterized by the bare mass m^, is given by 



= [i7-P+m,]e/+ J s,{e)g^K{p-e)':fj^ (2) 

(in the Euclidean space), while the bound-state vertex 
= Tgq' of the meson Ai composed of the quark q 
and antiquark q' , is 



(3) 
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In the above equations, integrations are over loop mo- 
menta, and Z)°^ (p) is an effective giuon pro pagator which 
at the present stage of DS studies p^. llOl. [20l . [2l1 | must 
be modeled in the non-perturbative QCD regime for low 
momenta, p^ %\ GeV^. 

In DS approach, all light pscudoscalar mesons 
(tt^'^, 77) manifest themselves both as quark- 

antiquark (gg) bound states and (almost-) Goldstone 
bosons of DChSB, which also generates the condensates 
and the weak decay constants of the right magnitude, 
rather independent of the small current quark masses 
fhq, which can even be vanishing, as fhq = corre- 
sponds to the chiral limit. For reviews, see, e.g., Refs. 
[il, [H, H Ref. [H] also reviews the studies of 
QCD DS equations at T > 0, started in 

The restoration of the chiral symmetry should take 
place at the temperature Tch, which is expected to be 
close or maybe even equal to the critical temperature 
Tc (e.g., see Refs. [Hill [11]). Most scenarios expect 
the restoration of the Ua{^) symmetry as the topolog- 
ical susceptibility melts with temperature. Such a sit- 
uation calls for a good understanding especially of the 
light pseudoscalar nonet, where pions and kaons to an 
excellent approximation, and 77 and 77' to a good approxi- 
mation, are eigenstates of flavor SU(3). However, for the 
neutral states with hidden strangeness, this will change 
if the gluon anomaly contribution melts with tempera- 
ture. The most interesting effects were predicted for 77 
and 77'. Their production rates were predicted to be en- 
hanced [H, [23, because their masses were expected to 
fall with T as the restoration of the Ua{^) symmetry 
takes place. However, for understanding such signatures, 
pseudoscalar meson masses (especially m,^ and tti^/) at 
T > must be understood better. As a contribution, 
the present paper extends to flnite temperatures the suc- 
cessful bound-state, DS approach of Refs. [1^, H^, [s^, [3l| 
to the 77-77' complex. 

The paper is organized as follows. In the next section 
we introduce briefly the chosen dynamical model, but 
also explain generally how in the DS approach one can 
construct the mass eigenstates in the 77-77' complex, be it 
at T = or T > 0. In the third section, the topologi- 
cal susceptibility is related to the anomalous part of the 
77 and 77' masses, and its temperature dependence fltted. 



Sec. llVl cxplains why various relationships between chiral 
and Uyi(l) symmetries would lead to different tempera- 
ture evolutions of the 77-7;' complex, and Sec. IVl gives the 
T-dependence of pseudoscalar masses for various scenar- 
ios. We summarize in the last section. 



II. PSEUDOSCALAR MESONS AT T > 

A. Model and its results in non-anomalous sector 

The non-Abclian ("gluon") QCD axial anomaly is es- 
sential for the 7^-77' complex. Refs. [H, [2^, [sO, HH show 
how to incorporate the effects of the anomaly into the 
DS-approach and achieve a successful description of 77 
and 77' (at T = 0). They also show [H, [H, [sO, [lH that 
this works if one first achieves what is needed for the 
no77-anomalous aspect of the 7;-7y' complex, namely the 
good description of pions and kaons. While the DS ap- 
proach in RLA tackles this efficiently at T = 0, it is 
technically quite difficult to extend solving Eqs. (Il|)-(l3|) 
to non- vanishing temperatures. We thus adopt a simple 
model for the strong dynamics, the separable interaction 



^iiuD{p^,q^,p- q) , 



(4) 



D{p\l\p-l) = DMp')M^^) + DiTi{p'){p-l)T^{e) .(5) 

This is an interaction with the two strength parameters 
Z?o, Di, and simply modeled form factors 



W) = 



.Fo(p2)=exp(V/A^), 
1 + cxp(-p2/A2) 



l + exp((p2_p^)/A^)- 



(6) 
(7) 



also used in the calculations at T > in Refs. |33l.l34l.l35l|. 
where fitting the properties of non-strange {NS) mesons, 
containing just u- and d-quarks, fixed the parameters to 

Ao== 758MeV, Ai = 961MeV, po = 600MeV, (8) 
Z?oA2 = 219, DiA;^ = 40, (9) 

which we also adopt here without any further re-fitting. 
As in Refs. [s^. Issl. Is^. Isst . the Matsubara formalism is 
used for calculations at T > 0. The usage of the above 
separable model interaction simplifies greatly DS equa- 
tions and calculations at finite T, while yielding equiv- 
alent results on a given level of truncation [33, |36|, and 
a similar quality of results for pions [s^ . l33l | and now 



also for kaons 
T = as in Refs 



(and fictitious ss pseudoscalars) at 
H, m, [33, HH , which employ more 
realistic interactions while modelin g Q CD. Most of the 
presently pertinent results of Ref. [3J| of the separable 
model at T = are summarized in Table HI In addi- 
tion, important are also the chiral-limit value of the qq 
condensate, {qq)^ ~ (— 217McV)"^ and the chiral-limit 
value of the pion decay constant, /° = 89McV. They 
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p 


Mp 




fp 


rexp 
Jp 


TV 


0.140 


0.1396 


0.092 


0.0924 ± 0.0003 


K 


0.495 


0.4937 


0.110 


0.1130 ±0.0010 


ss 


0.685 




0.119 





TABLE I: Results (in GeV) on the non-anomalous pseu- 
doscalar mesons for T — 0: masses Mp and decay con- 
stants fp of the pseudoscalar qq' bound states P = n,K 
and unphysical ss. The corresponding constituent masses are 
m„(0) = 0.398 GeV and m^O) = 0.672 GeV, the chiral- 
limit condensate is —{qq)o = (0.217)^ GeV^. All results are 
obtained for the bare quark masses mu,d = 5.5 MeV and 
nis = 115 MeV and the effective interaction parameter val- 
ues ((Sjl and ((9|, fixed in our Refs. [s^, H^. These masses 
and decay constants are the input for the description of the 
r^-rj' complex. [Later, in Eq. (|16|l. we will name the unphys- 
ical ss pseudoscalar meson rjs, but note that the mass M^g 
(|21|) . introduced in the NS-S mass matrix (I17|l . includes the 
contribution from the gluon anomaly, whereas Msa does not.] 



agree well with the values appropriate for QCD in the 
chiral limit [13, ■ The important aspect of the good 
chiral behavior is that the pseudoscalar masses obey the 
Gell-Mann-Oakes-Renner type relation 

M^^, = const (m, + mg>) , (10) 

as in all consistently formulated DS approaches, where 
Eq. (jlOp is an excellent approximation even for real- 
istically heavy strange quarks, so that the mass of the 
unphysical ss pseudoscalar can be expressed as 



MplGeVJ 




0.5 1 1.5 2 

FIG. 1: The dependence of the pseudoscalar meson masses 
on the relative temperature T/T^, where the temperature 
scale (characterizing rj-r^' behavior) is, for illustrative pur- 
poses, temporarily chosen to be (unrealistic) — 2/3 Tch = 
0.667 Tch- The chiral restoration temperature Tch is marked 
by the thin dashed vertical line. The masses which do not re- 
ceive contributions from the gluon anomaly, M^(r), Mk{T) 
and Mss{T), are depicted by the long-dashed, short dashed 
and dash-dotted curves, respectively. They practically do not 
change with temperature till T ~ 0.95 Tch, after which they 
rise monotonically towards the thin diagonal line, which (in all 
figures displaying masses) represents twice the zeroth Matsub- 
ara frequency, 27rT. This is the limit to which meson masses 
should ultimately approach from below at still higher tem- 
peratures, where qq states should totally dissolve into a gas 
of weakly interacting quarks and antiquarks. The masses in 
the rj-ri' complex exhibit different behavior due to the gluon 
anomaly, as will be explained below: the lower solid curve is 
Mrj{T), and the upper solid curve is M^i{T). The lower and 
upper dotted curves are, respectively, the rfj^g and rig masses. 



Mi 



2Ml - Ml 



(11) 



[E.g., the results in Table |T] obey Eq. (HH) up to \%.] 

Solving m Eq. ^ yields the masses Af,r, Mk and 
Mss and the BS vertices Tqqi . The meson decay constants 
fqqi (e.g., fus = Ik) can then be calculated as 



fqq' Pfj. 



(2^ 

X 757m Sq{i+^) Tqq,ii;p) Sq,{e~^) .(12) 



These calculations were extended to T > 0. To save 
space, the temperature-dependent masses of the pseu- 
doscalars tt, K and the ss are not displayed in a separate 
figure, but together with the anomalous masses from the 
r]-r]' complex, in Fig. [1] [and in Figs. [5]l9] further be- 
low, showing same Mt^{T), Mk{T) and Mss{T), just for 
variously scaled relative temperatures.] 

Fig. Ogives the T-dcpendcnccs of the decay constants 
for the realistic explicit chiral symmetry breaking, but 
also /°(T), the chiral- limit version of fn{T)- Its be- 
havior precisely shows the value of the chiral symmetry 
restoration temperature Tchj because f^{T) vanishes at 
T = Tch, as does the chiral condensate {qq{T))o. We 
also confirm [331 the result of earlier DS studies at finite 




TIGeVJ 



FIG. 2: The temperature dependence of the pseudoscalar me- 
son decay constants, /-^{T), fK{T) and /sa(T). We observe 
the typical crossover behavior. Even for the purely nonstrange 
pseudoscalar, the pion, f-rriT) does not fall to zero for the 
realistically small but non-vanishing explicit breaking of chi- 
ral symmetry, although the decrease of fiviT) is very strong 
around T = Tch, where it would even vanish in the chiral 
limit (i.e., if m, — > 0), as shown by the dashed curve repre- 
senting the chiral- limit pion decay constant /^(T). Increased 
explicit breaking of chiral symmetry due to s-quarks lessens 
temperature dependence, as shown by fK(T) and fss{T). 
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T [39j that the scalar cr-meson becomes degenerate with 
the pion after this temperature Tch; m the chiral Umit its 
mass Mcr(r) in fact tends to zero as T Tch [similarly 
as /°(r)], as should be at the chiral restoration. 

Besides all these correct, very advantageous behav- 
iors, there is the disadvantage that in the model used 
in Refs. (s^ . Issj and in the present paper, we find the 
value Tch = 128 McV, which is significantly lower than 
the values (167-188 MeV) indicated by various lattice 
QCD studies [13] ■ This is a well-known drawback of the 
rank-2 separable DS model [s^l, but its low chiral sym- 
metry restoration temperature can be raised by coupling 
Polyakov loop variable to the quark sector [4l|. This 
approach is currently under consideration [i^l • More im- 
portantly, note that in the present paper we do not ad- 
dress quantitative predictions at some specific absolute 
temperature; rather, as will become clear in detail in 
Sec. llVi we are interested in the dependence on the rel- 
ative temperature, since we will study various scenarios 
for various ratios between the chiral restoration temper- 
ature Tch and the temperature scale characterizing 
the disappearance of the gluon anomaly. For this, the 
present model is adequate in spite of relatively low Tch- 

In summary, the non-anomalous physics is all modeled 
successfully [sj,!!!] for the present purposes, and is ready 
to serve as input for constructing r] and 77'. 

B. The anomalous sector: 77—77' complex at T > 

The construction of physical meson states in the 
present model starts as in Refs. [H, IH, HO, [3l|. Be 
it at T = or r > 0, solving of the BS equations yields 
quark-antiquark bound-state solutions and the eigenval- 
ues of the squared masses; e.g., M^^^ = M|., as the pseu- 
doscalar us state is simply the positive kaon . By 
their strangeness and/or charge, such qq' states (where 
q' ^ q) are protected from mixing with the hidden fla- 
vor qq states uu, dd, ss, for which BS equations yields the 
masses M^^, M^j, M^j. Although the mass matrix of the 
light hidden flavor pseudoscalars is already diagonal in 
the basis uu-dd-ss, 

M%A = diag[ M„\, M^^-, M% ] , (13) 

the states mm, dd, ss do not correspond to any physical 
particles - certainly not to the isospin eigenstates 77^,77 
and 77', the neutral flavorless pseudoscalars found exper- 
imentally at T 0. Although M'^^ = M'^^ = in the 
isospin limit, A/^^, M^g, M'^g do not automatically repre- 
sent any physical masses, at least not at T = 0, since 
M^j^ is only the nan- anomalous (NA) part of the com- 
plete mass matrix = Mfj M\. Its anomalous part 
M\, the only nonvanishing one in the chiral limit, in the 
basis uu-dd-ss reads 



Ml = /3 
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breaking 
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X 


x^ 



where the first matrix neglects the influence of the SU(3) 
flavor symmetry breaking, while the second is its mod- 
ification which takes this into account as follows. Eq. 

says that due to the gluon anomaly, there are tran- 
sitions between the states mm, dd and ss. However, the 
amplitudes for the transitions from, and into, light pseu- 
doscalar mm and dd pairs are expected to be different, 
namely larger, than those for the significantly more mas- 
sive ss. To allow for the effects of the breaking of SU(3) 
flavor symmetry, we can write {qq\MWq'q') = bqbq>, 
where bg = ^ for q = u,d and bq = X^ for q ~ s. 
The most widely used estimate of the flavor breaking 
is the ratio of NS and S pseudoscalar decay constants, 
X = Uhs-s m, [m, m, m. since we compute both 
/tt and /ss, our A" is a predicted quantity, and not a 
fitting parameter. Our present model result at T = 0, 
/ttZ/ss = 0.773, is very close to A'cxp ~ 0.779 extracted 
phenomenologically, i.e., from the mass matrix featuring 
experimental meson masses - see Refs. [29l. [3l|. 

In any case, the complete mass matrix is very far 
from being diagonal in this basis, since phenomenology 
at T indicates [sec /3fit below Eq. dM])] that w 0.3 
GeV2 » Ml^, (and /3 ^ M^,). 

The basis suggested for neutral (/a — 0) mesons by the 
flavor SU(3) quark model, and especially by almost ex- 
act isospin symmetry, is the octet-singlet basis with well- 
defined isospin quantum numbers: the isovector (/ = 1) 
pion TJp = (mm — dd) j \/2, and the isoscalar (/ — 0) etas, 
778 = (mm 4- dd ~ 2ss)l\f^, 770 — (uu + dd + ss)/\/3. In 
the case of the flavor- broken SU(3), this is the effective 
octet-singlet basis, in which the complete mass matrix 
= M^j^ -\- M\ is rather close to diagonal in the fa- 
miliar T = case. That is, at zero temperature and for 
realistic flavor breaking, 77 is close to 778 and 77' to 770, 
thanks to the presence of the gluon anomaly (and in the 
chiral limit, the only non- vanishing mass is that of 770, 
then equal to 3/?). In the DS approach, this is shown 
in detail (with all numerical values) in Ref. [3l[. The 
present case is very similar, regardless the different dy- 
namical model used. 

Therefore, the effective flavor-broken SU(3) states 778 
and 7/0 will be close to the respective physical mesons 77 
and 77' as long as the gluon anomaly plays the similar role 
as at r = 0. However, vanishing of the anomaly {(3 = 0) 
at sufficiently high T would bring about the situation 
where the no-anomaly limit of the neutral flavorless pseu- 
doscalar states is: tt*^ = mm, 77 = dd, r( = ss (where even 
r]' would be an almost- Goldstone boson). This, in prin- 
ciple, opens the possibility to study the maximal isospin 
violation at high T [45|, [4q| , but as the effects of the small 
difference between u and d quark masses are not impor- 
tant for the present considerations, we stick to the isospin 
limit throughout the present paper. The pertinent states 
thus remain, even in the high-T-limit, the combinations 
with well-defined isospin quantum numbers. As long as 
MM and dd are bound, ir^ is given by (mm — dd)/^/2 and 
is decoupled from mixing with the two isoscalar etas, for 
which the no-anomaly-limit states are given by the so- 
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called NS-S basis of the 7 = subspace: 

1 - 1 [2 

riNs = ^("^ + dd) = —7=V8 + \ -^Vo , (15) 



V2' 



Vs 



V3' 
1 

71 



(16) 



The 2x2 submatrix of A'P 
the masses in the 77-77' complex, in this basis reads 



Mfj^ + M\ pertinent to 



Ml + 2/3 V2/3X 



A/2 
A4 



(17) 



where the indicated diagonalization is achieved by the 
NS-S mixing relations 

77 = COS07/JVS — sin(/)?75 , 77' = sin 077^^5 + cos i/i 775 . (18) 

The NS-S mixing angle 0, related to the (equivalent) 
effective 7^8-770 mixing angle 6 as 9 = (j) ~ arctan-\/2 = 
(p — 54.74°, is given by 



tan 2(j) 



][,p _ ]\,p 

VS VMS 



2V2(3X 

VS VNS 



where, from the 77-77' mass matrix (|17p . 



Mi 



VS 



Mi 



2/3, 



M^, + PX'^M^ + P 



J ss 



The theoretical 77 and 77' mass eigenvalues are 
1 



Mi 



M\ 



2 Wv.s 
1 



M, 



VS 



where A, 



(19) 

(20) 
(21) 

(22) 
(23) 



present paper we avoid the Ua(1) problem by breaking 
nonet symmetry just at the level of the masses, adding 
by hand the anomaly contribution M\ to the calculated 
non-anomalous mass matrix Af^^ to form the total mass 
matrix A/^. From the traces in Eq. pT)) . 



Mi 



Ml 



Mi 



(24) 



where M.^, Mgg and X — Jti/ fss are already calculated 
quantities (see Table |T]), whereas Af^ and M^i remain to 
be determined. 

One may treat /3 as a fitting parameter and fix it by 
requiring A/2 +Af 2 = {M^),^j, + {Mf^,),^p = 1.2169 GeV2 
in Eq. yielding /3fit = 0.282 GcV2. 

However, one can avoid treating /3 as a new fitting pa- 
rameter, since it can be extracted from the lattice. Using 
Eq. PT|) in Eq. gives the first equality in 
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13{2 + X')= M.^ + A/2 - 2A/^ = - X 

Jtt 



(25) 



The second equality in Eq. is the Witten-Veneziano 
(WV) relation [13, [3 between the 77, 77' and kaon masses 
and X: the topological susceptibility of the pure Yang- 
Mills gauge theory (x = Xym)- 

The lattice results (xiatt) on x imply through Eq. ([25]) 



/3i 



1 



latt 



6 

^2 Xlatt • 



(26) 



2 + 

Taking for xiatt the central value of the weighted average 

X(T = 0) = (175.7 ± 1.5 MeV)'^ , (27) 

of the recent lattice results on the topological suscepti- 
bihty [ii,!!^,!!! at T = 0, gives Aatt = 0.260 GeV2. 

Table summarizes our T = results on 77 and 77'. 
The agreement with the experimental masses is excellent 
for both /3 = /3fit and f3 = Putt- In the latter case, there 
is not even one new parameter introduced in addition 



III. TOPOLOGICAL SUSCEPTIBILITY AND 77, 

77' 

In the non-anomalous pscudoscalar meson sector, we 
calculated A/^r, Mgg, f-K, fss and X at zero and non- 
vanishing temperatures [3J|. Except for the anomalous 
contribution /3, this is everything one needs for com- 
puting, at T = and T > 0, the mixing angle ((T9|) 
and our predictions for the masses in the 77-77' complex. 
But, the anomaly effects, including /3, obviously can- 
not be calculated in the present DS approach utilizing 
the ladder approximation. Fortunately, since the gluon 
anomaly is suppressed as 1/Nc in the expansion in the 
number of colors Nc, it was shown that considering the 
gluon anomaly effect only at the level of mass shifts and 
neglecting its effects on the bound-state solutions is a 
meaningful approximation [2^, [2^, [3l| . Thus, also in the 



Mrjj^g Mrjg Mr,o M^, 

Pat 919.5 763.7 798.0 574.3 944.2 42.51° 548.9 958.5 

Aatt 884.0 735.2 790.0 573.7 914.8 40.82° 543.1 932.5 
Exp. 547.7 957.8 



TABLE II: The ?7 and rj' masses and mixing angle at T = 0, 
and comparison with experiment, for the cases /3 = /3flt = 
0.282 GeV^ and (3 = Aatt = 0.260 GeV^ For the both 
cases, the input from the non-anomalous sector is the same: 
X = 0.773 and and Mss from Table HI All masses M,,, 
M^i, Mn„^ s 8 0' ^^"^ are in MeV. The purely anomalous 
mass is the mass of 77' in the chiral limit for all three fla- 
vors, where rj is massless together with the other octet pseu- 
doscalars. In the chiral limit, as well as the exact flavor SU(3) 
limit, 4> = arctan v'2, i.e., 9 = 0, so that rj' = rjo, and r] = rja. 
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to the parameters already fixed in the non- anomalous, tt 
and K sector. The case /3 = /3iatt is especially satisfying 
because the lattice results for the topological suscepti- 
bility X exist also for T > [s^, [s^ , and can be used 
to express (3{T > 0) and study the 77-77' complex also at 
nonvanishing temperatures, if we assume that the WV 
relation continues to hold also at T > 0. 



Various lattice results [53 on the topological suscepti- 
bility x(r) are fitted well by the following expression. 



X(T) = x(0) cxp 



T 



(28) 

depending on the relative temperature T/T^, where 
is the characteristic ("melting") temperature where xiT) 
starts decreasing appreciably. 

For K = 4.2, a = 0.2 and 6 = 0.76 (solid curve in Fig. 
[S]), Eq. jMl) fits weU the lattice data [13] on the pure 
Yang-Mills topological susceptibility, x(^) = x{T)ym, 
which is the one consistent to use in the WV relation. 
However, we will also be interested in fitting the shape 
of xiT)/xiO) for the SU(3) quenched QCD data [fit ob- 
tained with K = 5.1, a = 0.5 and 5 = k/4 in Eq. 
dashed curve in Fig. [3] , as well as that shape for the case 
of the four-flavor QCD lattice data [H^l (dotted curve in 
Fig. m obtained with k = 13, a = 0.5 and S = 2.4). 

One extreme behavior of the topological susceptibil- 
ity, qualitatively described by Eq. p8|) for low values 
of T^, is the most traditional scenario, due to Pisarski 
and Wilczek ^ and studied, e.g., by Ref. [13 . Now, 
however, it appears to be disfavored, as then xC^) nielts 
away with T rather quickly, since one supposes that xiT) 
is due to the instanton contribution which decreases ex- 
ponentially as the so-called "Pisarski- Yaffe suppression 
factor" (53- Then, T^. would be significantly below Tch , 
which may be judged too unrealistic to consider [56|. 
Still, we consider it for illustrative purposes, since in 
the present bound-state approach in conjunction with 
xiT) = x(^)ym, one can get falling masses and enhanced 
r] and rj' production (predicted by, e.g., Refs. [H, \2T\ ) 
only in this case, and for = Tch but with x(2^) falling 
much more sharply than x(2^)ym, as will be shown below. 

The opposite extreme possibility is xi^) « const; i.e., 
Ua(1) symmetry is not restored around critical temper- 
aturCj but at much higher, possibly infinite T, e.g., see 
Ref. [53]. This possibility is approached by Eq. ([^5]) for 
growing T^. 

At zero temperature, rj^s and rjs are almost as "mis- 
aligned" with r] and ?/ as they can be, as shown by our 
NS-S mixing angle close to w 45°. This must be so 
for any successful model, since it is mandated by phe- 
nomenology at T = 0, as seen, e.g., in Ref. [sij. But, if 
the topological susceptibility vanishes when T becomes 
large enough, (3-^0 causes </> ^ in the end, along with 
V VNS,V' ^ Vs and M^{T), M^- Mss{T). 

Depending on how the mass contribution of the gluon 
anomaly varies with growing T, various interesting ef- 
fects can occur before this limit is reached. 




FIG. 3: The relative topological susceptibility x(T)/x{0) 
given by Eq. (|28|l fitting the lattice data [H^] for the three 
cases: pure Yang-Mills (solid curve, the slowest melting rate), 
SU(3) quenched (dashed curve, intermediate melting rate), 
and four-flavor QCD (dotted curve, the fastest melting rate). 



IV. CHIRAL VS. Ua(1) SYMMETRY 
RESTORATION AT T > 

The temperature evolution of the i]-tj' mass matrix (fT7|) 
results from the T-dcpendence not only of Xi but also 
of U, fss, and Msg. As shown by Eq. (gH]), the 
ratio x{T) I f-K{TY is crucial for the T-depcndence of the 
anomalous mass contribution. 

Thus, contrary to the assumptions often made in the 
literature (e.g., Ref. [11]), the decrease of the topolog- 
ical susceptibility does not imply automatic restoration 
of Uyi(l) symmetry (e.g., see Ref. [I^)- We will see 
that the assumption that it does, is correct only if x 
falls sufficiently faster than does, so that x{T) / f^{T) 
falls towards zero as T grows. We should thus explore 
what happens for various possible relationships between 
Ty. and Tch, the respective temperatures of the onset of 
the fast decrease of xC^) and of f-n{T). 

Many papers addressing the relationship of restoration 
of the chiral vs. Uyi(l) symmetry, e.g., Refs. [27l[58l. [soj . 
consider only one value of the temperature characterizing 
the melting of xi^)- Those of them [Hl,!!^] using the lat- 
tice results of Ref. [13] on the T-dependence of x, do not 
employ the rather high susceptibility melting tempera- 
ture (260 MeV) of that reference. They adopt the view 
of Ref. [5^ that this temperature {T-^ in our notation) 
is the critical temperature Tc, but rescale the lattice re- 
suhs m for x{T) to a lower value, T^^T^^ 150 MeV. 
(This in the ballpark appropriate to the full QCD, but 
significantly lower than the pure Yang-Mills case, where 
X consistent with the WV relation is computed.) Rescal- 
ing of is also relevant for us, since we are interested 
in the dependence on relative temperature, i.e., we study 
various relationships between and Tch- Below, we ex- 
hibit the case = Tch as well as the cases when is 
roughly 15% and 30% below and above Tch- 

Recall that the present DS model gave us Tch = 128 
MeV for the chiral restoration temperature. This is lower 
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than lattice results, since lattice finds that T( 
coincide with the critical temperature 



seems to 
2^,[2i] and 

presently estimates be between 167 and 188 MeV 

for 2+1 flavor QCD [40|. Nevertheless, our results have a 
more generic meaning than working in this simple model 
would indicate at the first glance; using a more realistic 
interaction in DS approach, such as in Ref. (39| . gives 
Tch above 150 McV, but the results obtained there oth- 
erwise look similar to the corresponding results here, just 
rescaled to the higher Tch- Notably, Ref. [3§| also found 
that /^(T) and M^(T), vanishing at T = Tch in the 
chiral limit, do not vanish any more when the explicit 
breaking of chiral symmetry is introduced, but exhibit 
a crossover, just as we find in the present model. Also, 
the pseudoscalar meson mass exhibits very little varia- 
tion below Tch just as the tt, K and ss masses here. 
The main change is thus pushing the fall-off behavior 
of /u(T'), which happens around Tch, to higher tempera- 
tures. The temperature-induced changes in the rj-r]' com- 
plex are thus given essentially by 2/3(T) ~ 4x{T)/f^{T) 
and 0{T)fl{T)/fUT) ~ 2x{T)/fUT) [see Eqs. m 
and (HI])]. Fig. [H shows that unlike f-^{T), fss{T) does 
not decrease by an order of magnitude across the exam- 
ined T-interval. Thus, the question what decreases ear- 
lier, x(r) or f^{T), and the related relationship between 
and Tch, is what causes important, qualitative differ- 
ences in T-evolution scenarios, and not what we happen 
to chose for our concrete dynamical model in the DS ap- 
proach. 

Let us here again note the advantage of the DS ap- 
proach for analyzing the hot QCD matter in which qq 
bound states and/or resonances persist beyond the criti- 
cal temperature marking not a real phase transition, but 
a relatively smooth crossover [lo| . The DS approach en- 
ables the meaningful computation of quantities charac- 
terizing such bound states or resonances. For example, 
note that Fukushima et al. [s^, who were the first to 
study the interplay of chiral symmetry restoration and 
melting of x in WV relation ([25]) , uses the Nambu-Jona- 
Lasinio (NJL) model (following, e.g., the pioneering pa- 
pers [1,0). However, the lack of confinement in the NJL 
model also leads to the disadvantage, important in the 
present context, that rj' is not bound at all, not even at 
T = 0. It decays into quarks and antiquarks due to its 
large mass, so that M^' in the NJL model is not a well 
defined quantity [s^. On the other hand, the present DS 
approach does not have such problems. 

Another advantage of the DS approach is the behavior 
of /7r(r). This approach has the correct chiral behavior 
also for finite T, so that fTr{T) falls markedly at T = Tch, 
but since the DS approach naturally incorporates also 
the effects of realistic explicit chiral symmetry breaking, 
/tt (T) stays non-zero [l^, HI] , albeit small, well beyond 
T = Tch- This is in contrast to the chiral-Hmit /°(T), 
where /"(T > Tch) = precludes its usage in WV rela- 
tion for T > Tch (e.g., Ref. [H or the studies [13, [H,!!! 
employing chiral Lagrangians). 



V. RESULTS ON 77-77' COMPLEX AT T > 
A. Yang-Mills topological susceptibility 

First, we explore the temperature dependence of the 
77-77' complex following from the usage of the pure Yang- 
Mills topological susceptibihty, x(^) = x{T)ym since it 
is the one consistent with the WV relation. One should 
keep in mind that the pure Yang-Mills case has the slow- 
est T-dependence, i.e., a significantly slower fall-off than 
cases when quarks are present (see Fig. [3]), which will 
motivate considering other cases in the subsection IV Bl 

For illustrative purposes, let us begin with an unrcal- 
istically low x-melting temperature, say T-^ = 2/3 Tch. 
This is a representative case when the topological sus- 
ceptibility x(2^) becomes very small well before Tch, i.e., 
well before f-n{T) can appreciably decrease. The tem- 
perature dependence of 77 and 77' masses for this case is 
shown in Fig. [TJ They remain practically as at T = 0, 
until T « O.QT;^, after which the giuon anomaly contri- 
butions start melting, and both AI^{T) and Mrji{T) start 
falling. After T « 1.2 7^ = 0.8 Tch, we have practically 
A/,,(T) = M,,^^(r) and M^/(T) = M^^(r) as 77 becomes 
pure rjffs and 77' becomes pure 77s, and even pure non- 
anomalous ss after T « 1.3 T^. After T = Tch, all masses 
rise, and 77 soon becomes degenerate with the pion. 

This behavior is reflected in the mixing angle 0(r) 
which for the choice = 2/3 Tch just follows the topo- 



logical susceptibility ([28]). falling monotonically from its 
T — value to zero - see Fig. (4] 

An important, potentially experimentally recognizable 
feature in Fig. [1] is the decrease of the masses of 77', 
and especially of 77. It starts already well before T = 
Tch, where the 77 mass is less than a half of its T = 
value. The decrease of the masses should result in an 
increase of their (especially 77's) multiplicities. This case, 
with significantly lower than Tch, thus illustrates well 




T[GeV] 



FIG. 4: The temperature dependence of the NS-S mixing an- 
gle (j) for = 2/3 Tch (dash-dotted curve), = 0.758 Tch 
(dashed curve, the case not displayed in other figures), = 
O.SSGTch (dotted curve), and — Tch (solid curve). In 
all the cases, the pure Yang-Mills topological susceptibility is 
used in the WV relation (gSj). 
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FIG. 5: The relative temperature dependence, on T/T^^, of 
the pseudoscalar meson masses for Ty^ — 0.836 Tch, i.e, Tch = 
1.20 marked by the thin dashed vertical line. The long- 
dashed, short dashed and dash-dotted curves again represent, 
respectively, M^(T), Mk{T) and Mss{T). The lower and 
upper solid curves are again M^{T) and M^i{T), respectively, 
and the lower and upper (except on a small segment) dotted 
curves are M^^g{T) and M^^{T). 



the T-dependences of the masses of pseudoscalars in the 
Pisarski-Wilczck scenario (MIHIIJI- 

Let us see how the behavior changes when is in- 
creased by, e.g., some 25%, to ~ 0.836 Tch- Unlike 
in Fig. [U T = Tch docs not mark the onset of the ap- 
proach to the asymptotic situation [Mn'{T) Mss{T), 
M^{T) M^{T)]. Now, in Fig. [3 T = Tch marks the 
strong increase of the rj' mass, which almost doubles by 
T/Ty. w 1.4, where it reaches its maximum and starts 
fahing till its second, very narrow anti-crossing with the 
rj mass, which just grows tamely and (almost) monotoni- 
cally. [The crossings of 777^5 and 775 correspond to the anti- 
crossings of the 77 and 77' masses ([22]) and ([23]) .] The in- 
crease of Ty also pushed the onset of the approach to the 
asymptotic situation beyond this second anti-crossing. 

This all happens because now the decrease of the pion 
decay constant /^(r) starts sufficiently before the de- 
crease of the topological susceptibility, so that a large am- 
plification of (i{T) ~ x{T) / f-^{T) occurs, and over a sig- 
nificant temperature interval. In contrast, /3{T)X(T)'^ ~ 
x{T) / fss{T) is not so significantly enhanced. Eqs. pO]) . 
(PT|) and (pSj) show why M^^^ is larger than M^^ in the 
interval 1.1 ^ T < 1.75 and is the largest contrib- 
utor to Alrj'. In the T- interval where Mr,^g{T) is signifi- 
cantly larger than other masses and /3(T)X(T)^, expan- 
sion of A,pj' in Eq. ([22|l shows that the 77 mass is approx- 
imately the non-anomalous ss mass Mss{T), with the 
gluon-anomaly contributions suppressed by l/AiriNsi"^)'- 



Mi 



ss 2 



n 2 



VMS 



M, 



Vns 



+ 



(29) 



In this temperature interval, the heavier particle, 77', is 
predominantly non-strange, rjpfs, while the lighter one, 
77, is mostly strange, 775. Eq. ([TO)) shows that as AI^^^ 
grows, so docs 0(T). It surpasses 45° when Airing sur- 



passes M^^ - compare Figs. |4] and [5l In fact, after 
T « 1.25 = Tch, it climbs close to 0(T) ~ 90°, as 
77' 77Ar5 and 77 ~ -775. After, as x(T) -> 0, A/,,„s(r) 
falls back towards M^(r), and (^(T) falls [through 45° 
where again AI^ihs — Mjj^, while A'/^'(r) and A/^(r) 
anti-cross] towards zero, and etas quickly reach their no- 
anomaly limit (jl5|) - (fTCl) : r/ = tins degenerate with the 
pion, and 77' = 775 with the non-anomalous mass Mss{T). 

The crucial role of the chiral restoration tempera- 
ture Tch is already obvious. We assume next that the 
topological-susceptibility-mclting temperature is still 
higher and equal to this temperature, i.e., = Tch. 
This case is illustrated in Fig. [6l It exhibits several in- 
teresting features. The most obvious one is such a strong 
increase of the 77' mass that soon after T = Tch, it ex- 
ceeds the scale of the figure. Of course, it does not indi- 
cate any divergence; after reaching the maximum value of 
around 5 GeV, M^y [T) falls down. It is just that the de- 
creasing part of the A/,j/(r) curve below 1.6 GeV and the 
close r;-7;' anti-crossing are pushed outside the examined 
temperature interval by widening of the interval where 
MriNs (^) is enhanced. This results in the mixing angle 
(j)(T) which in the displayed T-interval (T < 200 MeV), 
does not return down from the value of approximately 
90°; i.e., 77' remains almost pure 77JVS, and 77 almost pure 
775 in the displayed T-interval. The fall of M^/ and ap- 
proach to the no-anomaly asymptotic situation is thus 
postponed to higher temperatures now that we increased 
Ty^, the 'melting temperature' of the (anyway compar- 
atively slowly decreasing) Yang-Mills x- On the other 
hand, since the decrease of f-K{T) is determined by Tch, 
increasing T-^ lowers the value of the relative temperature 
where M^^^. (T) begins its fast climb. 

All these features are readily understood on the ba- 
sis of the detailed explanation of the previous case = 
0.836 Tch, but they arc much more pronounced with re- 
spect to this case with just 16% lower T^. The described 
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FIG. 6: The relative temperature dependence, on T /T^, of 
the pseudoscalar meson masses for = Tch, marked by the 
thin dashed vertical line. The meaning of all curves, lines and 
other symbols is the same as in the Figures[T]and[5l The most 
conspicuous is the increase of AI^i{T), which become as high 
as 5 GeV around T = 1.3 T^. 
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FIG. 7: The T /T^ dependence of the meson masses for T-^ — 
1.17 Tch- The meaning of all symbols is again the same as 
in Figs. [1] \5\ and El The increase of M^r (T) is huge, as high 
as 10 GeV around T = 0.18 GeV, since the masses in the 
rj-Tj' complex are calculated (as in Figs. [T] [5] and [6] where 
Tx is, however, lower) with the pure Yang-Mills topological 
susceptibility [S^, which decreases comparatively slowly. 



trends continue for > Tch- We illustrate this through 
the choice = 1.17 Tch, which in our model corresponds 
to = 150 MeV. This value is equal to T-^ and the crit- 
ical temperature assumed, e.g., in the typical examples 
[27l . [ssl . [sol of studies paralleling ours. These assump- 
tions are obviously motivated by 150 MeV being roughly 
the value of the critical, transition temperature Tc found 
on the lattice for 3- flavor QCD 

The temperature behavior of the pseudoscalar meson 
masses for this case is shown in Fig. [71 and it is clear 
from the understanding gained on the previous choices 
of . The increase of the rj and 77' masses starts as early 
as T/T^ — 0.6, and, in the case of the rj' mass, it is even 
more pronounced, the maximal Af^/ (T) being some 10 
GeV. Its fall towards the ss mass (and the fall of the 
mixing angle from almost 90° to 0) is not seen in Fig. |7| 
since the approach to the asymptotic behavior happens 
at higher temperatures. 

Let us summarize our findings from the WV relation 
used consistently with the pure Yang-Mills topological 
susceptibility: only very low, unrealistic values of T^, be- 
low some 70% of Tch, would yield the drop in the p and 
T]' masses which would allow the enhancement [2^, of 
the T] and 77' production as in Pisarski-Wilczek scenario. 
For all higher values of we found mostly increasing 
r] and 77' masses allowing no such enhancement in their 
production. The increase of the 77 mass in Figs. |5l |6l 
and [71 is however moderate and relatively slow, so that 
in our qualitative treatment wc do not venture to make 
any statement about a possible suppression of the rj pro- 
duction either. In contrast to that, the enhancement of 
the r]' mass is so dramatic for > Tch , that one can ar- 
gue that the 77' production must be ever more suppressed 
for these values of - of course, provided that the 77' 
mass increase observed in the present approach is at all 
genuine, and not an artefact. 



B. Discussing caveats and probing other 
possibilities 

While our bound-state approach, which includes both 
DChSB and realistic explicit ChSB, makes it possible to 
use WV relation, and the pion decay constant in it, for 
the finite temperatures across the restorations of Ua{^) 
and chiral symmetries, there is no guarantee that the WV 
relation does not ultimately fail at some T. Of course, it 
is not likely that the WV relation, after being a good ap- 
proximation at lower T, should fail almost immediately 
after f^r starts falling with T. Thus, the behaviors of 77 
and 7/ exhibited in the previous subsection should be con- 
sidered at least as a serious qualitative indication that, 
e.g., some 77' mass enhancement (and associated suppres- 
sion of 77' production) should exist at least over some T- 
interval, even though the WV relation may break down 
after certain T. At this point, nevertheless, we must dis- 
cuss and take into account the caveat that such a drastic 
77' mass increase (for > Tch) suggest that it is in fact 
not a genuine effect, but at least to a large extent an 
artefact of our approximations and assumptions. 

One can think of the large Nc approximation, in which 
the WV relation (with x — Xym) was derived, as a pos- 
sible cause of unreliability of the WV relation after some 
T. The avenue towards improvement and better insight 
can then be sought in the recent work by Shore [61] . [62 | , 
which contains what amounts to the generalization of 
the WV relation vahd to all orders in 1/iVc- Instead of 
the topological susceptibility x, Shore's equations employ 
the full QCD topological charge parameter A (but since 
it is at present not known. Shore himself (6l] . had to 
approximate it by its lowest-order approximation in the 
1/Nc expansion, namely by the pure Yang-Mills topolog- 
ical susceptibility xym)- Also, besides /tt, the kaon decay 
constant and four different 77 and 77' decay constants ap- 
pear in the pertinent Shore's equations (2.14)-(2.16) [6]| . 
Very illustrative is also Eq. (2.19) in Ref. [6y, which af- 
ter the 1/Nc expansion gives back the WV relation ([25]) 
as the lowest-order approximation. So, instead of the 
WV relation. Shore's generalization might be used in our 
bound-state approach. Indeed, this research has been 
initiated and is in progress [42|, and there are already 
some preliminary results (notably at = Tch, but the 
presentation of this material is beyond the scope of the 
present paper. Nevertheless, already on the basis of the 
insights from the present work we can and should make 
some comments on the indicative fact that the results 
from Shore's generalization at T > arc rather close [1^ 
to our present results from the WV relation provided the 
same topological susceptibility % is used instead of the 
full QCD topological charge parameter A. 

It is instructive to compare the WV relation ((25|l with 
Shore's [H Eq. (2.19). The refinement brought by the 
presence of the five additional decay constants in Shore's 
generalization (instead of just as in the WV relation) 
immediately suggests that the rj' mass increase would not 
be as drastic as the ones obtained from the original WV 
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FIG. 8: The relative temperature dependence, on T /T^, 
of the pseudoscalar meson masses for the SU(3) quenched 
topological susceptibility from the lattice [13] • From top 
to bottom, the diagrams correspond to = 0.667 Tch, 
= 0.836rch, Tx = Tch, and = l.lTTch. That means 
the diagrams in this composite figure correspond to the in- 
dividual Figs, [l] m in and [7] (employing the pure Yang-Mills 
topological susceptibility). The description of the curves and 
symbols is the same as in these figures. In the third and fourth 
diagram, M^i exceeds the scale, and the respective maxima 
are roughly 2.8 GeV and 7.7 GeV. What looks like a crossing 
of Mrf and M^i in the third diagram, is really a very close 
anti-crossing, as clearly seen in the third diagram of Fig. |9] 



relation for T^/Tch > 1. Namely, these additional decay 
constants are all affected by the strange quarks; the way 
fss contributes to the ry and 77' decay constants can be 
seen in, e.g., Ref. [i^ or Appendix of Ref. [2^. Thus, 
they diminish with T slower than as shown in Fig. 
[5] for fx and fag- In keeping with that, our preliminary 
results indeed show that the 77' mass increase is not as 
pronounced as in the present paper. However, the re- 
duction is just about some 20%, and, as abeady said, 
the results are on the whole rather similar. The refine- 
ment through the decay constants (other than /^r) thus 
improves, but does not really cure the suspicious yy'-mass 
behavior. This indicates that the probable main problem 
is the usage of the topological susceptibility x = Xym 
instead of the presently unknown full QCD topological 
charge parameter A. While the usage of the pure Yang- 
Mills topological susceptibility xym often turns out to 
be a reasonable approximation at T = (e.g., see Refs. 
[H, m, S im, .43, 61, 6^]), its 'slow' temperature de- 
pendence contributes crucially to the perceived artefact 
of the present approach - the blow-up of the ry' mass. 
This motivates us to use the related topological suscep- 
tibilities obtained on the lattice [s^ in the presence of 
quarks, namely the SU(3) quenched and the four-flavor 
QCD case, as Fig. [3] shows that they 'melt' faster. Thus, 
they hopefully provide a better approximation to the T- 
dependence of the full QCD topological charge parame- 
ter, since we expect that, as usual, this quantity of the full 
QCD will be less resistant to the temperature-induced 
changes than the pure-glue, quarkless xym- 

In Fig. [8] we see the results for the SU(3) quenched 
QCD data, and in Fig. [9l the results for the four-flavor 
QCD data [as fitted by, respectively, the dashed and the 
dotted curve in Fig. [3l resulting from Eq. ([28])]. 

In the both figures, the same four values of the relative 
susceptibility melting temperatures T^/Tch are used as 
in the figures in the preceeding subsection. Thus, in both 
Figs. [8]and[9l the four successive graphs starting from 
the top, correspond to the four respective Figs. [U O [H 
and[7l which gave the T-dependence for x = Xym- It was 
described in detail in the previous subsection, and on the 
basis of this we can easily understand the behavior of the 
masses in the new figures [8| and [9l What happens is that 
for the SU(3) quenched QCD, and especially for the four- 
flavor QCD, the topological susceptibifity melts faster 
than in the pure Yang-Mills case, so that the patterns 
seen in the previous subsection for the pure Yang-Mills 
case are (roughly) reproduced at somewhat higher values 
of for the SU(3) quenched case, and even higher values 
of for the four-fiavor case. Thus, the Pisarski-Wilczek 
scenario is realized in the upper two diagrams of Fig. 
[9] illustrating the four-flavor case, and just in the flrst 
diagram of Fig. [8l but in any case we observe Af,,(T) 
become degenerate with Mt^{T) fast, before the chiral 
restoration temperature Tch, while the behavior observed 
previously in Fig. [T] would now be still seen well above 
T^/Tch = 2/3. The second diagram from the top of Fig. 
[S] illustrates nicely how the bump in .A/.^„g starts growing 
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but then fails to develop since now x(^) nielts too fast 
at T^/TcY^ = 0.836, in contrast to the mass bump in Fig. 
[51 That bump, namely the strong enhancement of the rj' 
thermal mass, is now for the both susceptibilities, i.e., in 
both Figs. [8]and[9l encountered in their third diagram, 
where T^/Tch = 1- For the four- flavor QCD, with its 
faster melting x{T), the enhancement is only some 60%, 
but for the SU(3) quenched QCD [with an intermediate 
rate of the x(T) melting], the enhancement is drastic so 
that and M^„g mass bumps exceed the displayed 
scale of the figure, the maximum being roughly Af,,/ « 
2.8 GeV. For the four-flavor QCD, such a mass bump 
happens in the last, "highest-T^/Tch" diagram of Fig. [9l 
The last diagram in the SU(3) quenched case. Fig. [51 of 
course exhibits even stronger ilf^/ enhancement, although 
weaker than in Figs. [Sland[71in the pure Yang-Mills case. 

In summary, we tested the sensitivity of our results on 
variations of topological susceptibility, and found that 
they are sensitive to the detailed temperature behavior 
of x(2^)- The changes are readily understood as the 
consequences of the faster-melting x's- The relative- 
temperature behaviors of the 77 mass, both in this and the 
previous subsection, is consistent with the M^(T) behav- 
iors recently found in the NJL model by Costa et al. [6^ 
for various melting rates of x- The enhancement of the ry' 
mass remained the most conspicuous feature for the real- 
istic cases > Tch, but it was reduced somewhat by the 
faster-melting susceptibilities. Although most often still 
suspiciously large, the M^/ (T) enhancement looks reason- 
able for the case of the four-flavor QCD with = Tch 
(the third diagram in Fig. [9]), where one should keep in 
mind that it would be further reduced in Shore's gener- 
alization. 



VI. SUMMARY AND CONCLUSIONS 

The separable model of Ref. [32| provides the good 
description of the nonstrange light mesons at T = and 
T > 0. We extended it to the strange sector: the kaon 
phenomenology was successfully obtained, along with the 
results on the fictitious ss pseudoscalar meson, needed 
for the subsequent calculations of the rj-rj' complex. The 
rj-r]' phenomenology at T = was then also successfully 
reproduced, now in this specific dynamical model [s^l, by 
following Refs. [H, [2^, [s^, [3l|. Thereby the stage was 
set for the T > extension of the DS approach to rj and 
77' [H, [2^, [33, HH , which is the main goal of the present 
paper. 

The WV relation ([25]) enables the DS approach f28l.[29l. 
l30l . Isij l to determine the masses in the 77-77' complex in a 
parameter-free way, from the calculated non-anomalous 
pseudoscalar qq masses, the calculated pion decay con- 
stant /tt (resulting from our qq bound-state solutions) 
and the lattice results on the topological susceptibility 
X, be it at T = or T > 0. 

The results for the thermal evolution of the masses 
in the rj-rj' complex thus depend on the choice of the 




FIG. 9: The relative temperature dependence, on T/Tx, 
of the pseudoscalar meson masses. The masses in the r^-rj' 
complex are calculated employing the topological suscepti- 
bility for the four-flavor QCD taken from the lattice. From 
top to bottom, the diagrams correspond to — 0.667 Tch, 
= 0.836 Tch, Tx = ^Ch, and = l.lTTch. That is, the 
diagrams in this composite figure correspond to the paral- 
lel diagrams in the previous figure [obtained with the SU(3) 
quenched x{T)], as well respective as the individual Figs. [T] 
[5l[6]and[7] (obtained with the pure Yang-Mills topological sus- 
ceptibility) . The description of the curves and symbols is the 
same as in these figures. In the fourth diagram, M^/ exceeds 
the scale, and its maximum is 6.4 GeV. 
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topological susceptibility x(T); in addition, the results 
depend even more on the relation between the chiral 
restoration temperature Tch and T^, the temperature of 
melting of the topological susceptibility. Let us recall the 
oldest considered scenario, namely the Pisarski-Wilczek 
one, where the expected drop of the 77 and 77' masses 
as x(T) drops with T, led to the expectations of the 
enhancement of ry and 77' multiplicities with growing T 
[26l . \2T\ . The (relative) T-dependences of the masses con- 
sistent with this scenario are realized in Fig. [U the first 
diagram of Fig. [8] and the first two diagrams of Fig. [9l 
However, this requires a rather low T^, noticeably below 
Tch- This is unrealistic, as the values of lower than 
Tch are excluded by the lattice; e.g., AUes et al. have (in 
their notation, T^, T^) « 260 MeV [H, [6i|. They 
identify this with Tch, but for the value appropriate 
to the pure Yang-Mills case, which is higher than Tch in 
the presence of quarks. Also, there are the analogous, 
newer lattice results [131, with even higher « 300 
MeV. In the both cases, is the melting temperature 
of the topological susceptibility of the pure gauge, Yang- 
Mills theory. It is thus natural that it is above the char- 
acteristic temperatures of the full QCD, which are low- 
ered (with respect to Yang-Mills) by the presence of the 
quark degrees of freedom. This may cause (in analogy 
with Tch) that from the high values characterizing the 
pure Yang-Mills case, the physically relevant (e.g., 
appropriate to the full QCD topological charge parame- 
ter) gets lowered to Tc ~ Tch of the full QCD - but not 
to still lower values (T^ < Tch) generally required to re- 
alize the Pisarski-Wilczek scenario in our model. Thus, 
our model indicates that this scenario and the hot QCD 
matter signal (increased 77 and 77' multiplicity) associated 
with it [29 . W?\ is excluded - apart from a possible partial 
exception, the third diagram in Fig. [51 This is the only 
> Tch case where we observe that M,^ (T) suffers a sig- 
nificant fall. After T = 1.1 Tch, it drops to roughly half 
of its T = value, i.e., to Mt^{T) which is not yet exces- 
sively thermally enhanced. Being preceded by a (modest) 
rise, this fall of M^(T) maybe cannot provide an easily 
detectable increase of the 77 production, but this is closest 
our DS approach gets to the Pisarski-Wilczek scenario for 
an acceptable value of T^/Tch- 

Having mentioned the full QCD topological charge pa- 
rameter appearing in Shore's generalization of the WV 
relation, let us add that the present paper and issues risen 



therein, obviously provide additional motivation for lat- 
tice calculations to determine this quantity and its melt- 
ing temperature T^. We expect that this may turn 
out to be equal or close to the other characteristic tem- 
peratures of the full QCD, in analogy with the result 
Tc « Tch; e.g., see Refs. fUlHlli. 

This lowering of to Tch would help controlling the 
77' mass enhancements appearing for ^ Tch, which 
become more and more dramatic as T^/Tch grows. Ex- 
cept in the third diagram of Fig. [HI these enhancements 
are so drastic, that they must to a large extent be an 
artefact of our approximations and assumptions, and one 
cannot rely on them quantitatively. On the other hand, 
their persistent appearances are at least a qualitative in- 
dication that probably M^/(r > Tch) > Mrj'{0) in some 
T- interval around Tch. Thus, even though we cannot 
make quantitative predictions, on the basis of our results 
we expect a suppression of the rj' yield after the onset of 
the chiral symmetry restoration in the hot QCD medium. 
We thus propose that the measurements of the rf mul- 
tiplicity should be undertaken at RHIC and LHC. Even 
if the 77' suppression would not be seen, it would still be 
an interesting result because it would falsify the T > 
extension of the so far successful bound-state approach 
to 77 and 77' [2i,[2§,[30,[3li. 

In contrast to 77', the mass of 77 does not exhibit a 
marked rise in any figure (except the gradual approach 
to the limit of two lowest Matsubara frequencies). In 
fact, beyond the chiral restoration temperature, it rises 
somewhat more slowly than the pion mass. On the ba- 
sis of these masses, one cannot expect any suppression 
of the relative tj/tt'^ multiplicity. This is in agreement 
with the recent experimental results of PHENIX collab- 
oration, which recently found a common suppression pat- 
tern of 77 and tt" mesons at high transverse momentum 
in Au-|-Au collisions [65| . 
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